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Abstract 

We analyze the effects of flavor splitting from staggered fermion lattice simulations of a low 
energy effective theory for graphene. Both the unimproved action and the tadpole improved action 
with a Naik term show significant flavor symmetry breaking in the spectrum of the Dirac operator. 
Note that this is true even in the vicinity of the second order phase transition point where it has 
been argued that the flavor symmetry breaking should be small due to the continuum limit being 
approached. We show that at weaker couplings the flavor splitting is drastically reduced by stout 
link smearing, while this mechanism is ineffective at the stronger couplings relevant to suspended 
graphene. We also measure the average plaquette and describe how it calls for a reinterpretation 
of previous lattice Monte Carlo simulation results, due to tadpole improvement. After taking into 
account these effects, we conclude that previous lattice simulations are possibly indicative of an 
insulating phase, although the effective number of light flavors could be effectively less than two 
due to the flavor splitting effects. If that is true, then simulations with truly chiral fermions (such 
as overlap fermions) are needed in order to settle the question. 
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I. INTRODUCTION 



Recently, a number of lattice Monte Carlo simulations of graphene and graphene-like sys- 
tems have appeared [l-4]. Refs. 1-3] study the effective theory of Nf flavors of massless four 
component Dirac fermions, constrained to 2+1 dimensions, subject to an instantaneous 3+1 

nn n 

dimensional Coulomb interaction [5|, |6[. In Ref. |4J, a 2+1 dimensional Thirring-like model 
is investigated. This is related to graphene-like systems through a large Nf or strong cou- 
pling equivalence in the dispersion relation for the auxiliary boson versus photon. Graphene 
has Nf = 2, but studying other Nf is of interest in order to understand the phases of 
such theories more generally, and because the large Nf limit is under theoretical control 
6|. Other interestingstudies coming from the effective field theory perspective have also 
recently appeared (JS]- 

In this article we address the flavor symmetry breaking that is introduced when staggered 
fermions are used in the lattice formulation. We also discuss the effect of photon tadpoles 
that come from lattice field theory. We will show that both features play an important 
role in the interpretation of lattice results. We explore various improvements to the lattice 
formulation. One is adding a Naik term to the action, which reduces discretization errors 
from 0(a) to 0(a 2 ), where a is the lattice spacing. Another is tadpole improvement, which 
removes ultraviolet divergent renormalizations associated with the lattice link operators. A 
final improvement that we consider is stout link smearing, which we find restores flavor sym- 
metries at weak couplings but not at the strong couplings relevant to suspended graphene. 
Importantly, we find that flavor symmetry breaking is significant in the vicinity of the sec- 
ond order phase transition point that occurs in the noncompact gauge formulation. Thus 
although it has been argued js] that the continuum limit should be approached at this point, 
and hence flavor symmetry violations (which are 0(a)) should be small in this regime, we 
have empirical results which contradict this expectation. Finally, we discuss how the flavor 
symmetry violations, revealed in split eigenvalues of the Dirac operator spectrum, perhaps 
imply that there are effectively less light flavors than two. Given the phase diagram that has 
been suggested by a number of studies in the Nf versus inverse coupling plane, this would 
imply that the critical coupling for Nf = 2 would occur at a somewhat stronger coupling 
than is found from staggered fermions. Only a simulation with truly chiral lattice fermions, 
such as overlap (Neuberger) fermions Jq], can conclusively answer the question of what is the 
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critical coupling for Nf = 2, since no systematic way of restoring the flavor symmetry has 
been found so far for the staggered fermion formulations at the stronger values of couplings. 

The outline of this paper is as follows. In Section [TT] we describe the action of the 
continuum effective theory that is supposed to describe the low energy limit of suspended 
graphene. We pay particular attention to redefinitions that are involved in going to the 
action in its simplest form, as these will be mirrored in redefinitions made in the lattice 
formulation. It will be shown that in the massless limit there is only one parameter in the 
theory, a coupling g which is strong in the case of suspended graphene. We also describe the 
£7(4) flavor symmetry of the effective theory, which is spontaneously broken to {7(2) x {7(2) 
by the formation of a chiral condensate, when the coupling g is sufficiently strong. In 
Section Hm we discretize the continuum action, formulating the lattice theory with staggered 
fermions. We show the redefinitions that isolate the one parameter of the lattice theory (in 
the massless limit), (3 = 1/g 2 . It is important here that we make redefinitions that maintain 
the unitarity of links; i.e., U(n) = exp[i6(n)], where 6{n) is a real lattice field representing the 
scalar potential associated with the instantaneous Coulomb interaction. Interestingly, this 
approach demands an anisotropic lattice with lattice spacing a t in the time direction and a s 
in the spatial directions, with the anisotropy parameter a s /a t set equal to the Fermi speed, 
a s /at = vp. Flavor symmetry violation of the unimproved staggered fermion formulation 
is discussed in Section HVl For 2+1-dimensional staggered fermions, 0(a s ,a t ) terms reduce 
the £/(4) flavor symmetry to {7(1) x {7(1) in the massless limit. We evaluate the spectrum 
of the unimproved Dirac operator on a large number of lattice field configurations that we 
have generated by Monte Carlo techniques. We show that at stronger values of the coupling 
g (equivalent to small values of /3), the flavor symmetry violation is severe. This is revealed 
by the lack of four-fold spectrum degeneracies that would be present if the {7(4) symmetry 
were respected. We find that this is even true near the second order phase transition point 
of the noncompact gauge formulation. 

Section |V] describes Naik fermion and tadpole improvements to the lattice formulation. 
We do not find any restoration of flavor degeneracy but do find significant reinterpretation of 
the bare lattice parameters in terms of those that are tadpole improved, at strong coupling. 
It will be seen that that has important implications for the phase diagram of the theory. 
We measure the average plaquette in dynamical simulations. We will show that for stronger 
couplings the resulting tadpole improvement of the theory has a large effect when relating 
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the simulation lattice coupling (5 to the coupling in the tadpole improved action, /3ti- The 
result is that for the noncompact gauge action the insulator /semi- metal transition occurs 
at a physical coupling that is significantly smaller than the g 2 of suspended graphene. The 
apparent absence of a spectral gap in the experimental results for suspended graphene near 
the Dirac K points [14j is in conflict with the lattice simulations, and we will not be able to 
provide an explanation for this discrepancy. 

The topic of stout link smearing, which is also a type of improvement, is discussed in 
Section ED We find that this is very effective at weak couplings, but that it is not useful for 
restoring flavor symmetry at the strong coupling relevant to either graphene or the second 
order phase transition point that occurs in the noncompact gauge action. We conclude in 
Section IVIII with a number of observations, summarizing our finding. 



II. CONTINUUM ACTION 
A. The effective coupling g 

The Euclidean spacetime action for the effective theory is given by 

S = dtd 2 X ^ {tpaloDflpa + flV F ^ i>alidii) a 
J a=l,2 i=l,2 

W^a) + | / dtd 3 X ( d M 2 (1) 

i=l 

Here 7$, i = 0,1,2, are Euclidean Dirac matrices satisfying the 50(3) Euclidean rotation 
group Clifford algebra {7*, 7j} = 25^. For instance we could choose 

( 

7* = , i = 0,1,2, (2) 

\-m J 

composed of Pauli matrices with <7o = (73. Also note that due to the nonrelativistic approx- 
imation, the covariant derivative only involves the scalar potential A 

D t = hd t - ieA (3) 

Next we make the redefinitions 

x = v F t, A = —^A' Q} D = d - iA' (4) 
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to obtain 

= f d 3 X (TpaJoDolpa + ^ ^ali^c 

a=l,2 i=l,2 

F ^ 8=1 

Recall that the Euclidean path integral that defines the theory has as its integrand 
exj)(—S/h). The rescalings have isolated the sole coupling constant in the theory, 

e 2 

9 2 = t = {c/v F )4na (6) 

where a is the fine structure constant. A final redefinition 

A' = gA , D = d - igA (7) 

makes it clear that g is the coupling constant in the photon-electron-electron vertex of this 
theory. 

2 

Perturbation theory would be valid in the limit where a g = |- = a c/vp <C 1, which 
is clearly not the case for graphene, where c/vp ~ 300. Given that the coupling is in fact 
strong, it is natural to appeal to lattice Monte Carlo methods, as has been done in the 
case of the nuclear strong interaction, quantum chromodynamics (QCD). It also becomes 
clear why one would like to be able to adjust vf experimentally, since the coupling of the 
theory determines the binding energy of any possible bound states that might form from the 



massless quasiparticles, analogous to hadrons in QCD. In fact, [l-|4| argue that the theory 
is quite similar to QCD in that when the coupling is strong enough one creates a nonzero 
"chiral" condensate (ipaipp) ^ 0, so that the theory is in a Mott insulator phase. (Properly 
speaking, chirality does not exist in 2+1 dimensions. It is, rather, a flavor symmetry that 
is being spontaneously broken in the 2+1 dimensional effective theory.) 

B. Symmetries 

The three-dimensional 50(3) ~ SU(2) rotation group acting on the spinors has genera- 



tors Sij = \oij ® 1 where 



°ij = -(V 2 )[7i ? lj] = tijk diag (cr fc , a k ) (8) 
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and the 1 factor in \(J{j <8 1 acts on the two dimensional flavor space. The action (j2J) has a 
£7(4) flavor symmetry, with 16 generators that commute with those of the rotation group, 

1(8)1, 1 <8 <Ti, 7475 <8 1, 7475 <8 cr* (9) 
74 <8l, 74 <8 (Ji, 75 (8 1, 75 <8 <Ji (10) 



where 7 4j5 are given by 





74 = ,75 = (11) 



when we choose the Dirac matrices (j2J) • A mass term m Yl a ipaipa reduces the symmetry 
to U (2) ® U(2) since the generators ( flOi) are broken. However, we still expect a four-fold 
degeneracy in the spectrum of the Dirac operator 

,2 

hvi 



M = l0 D + + TT~ (I 2 ) 



i=l,2 

because the 4 representation of Z7(4) decomposes to a (2, 2) representation of the subgroup 
SU{2)®SU{2). (In spin language, this is the (j u j 2 ) = (1/2, 1/2) representation of SU(2) 1 ® 
SU{2) 2 -) This is important in our considerations below because the Monte Carlo simulations 
are done at a nonzero mass, in order to avoid numerical difficulties (inversion of a poorly 
conditioned Dirac matrix). We will examine the spectrum of the Dirac operator on the 
lattice and compare to this four-fold degeneracy of the continuum theory with a mass term. 

The formation of a "crural" condensate ( / ip a ipp) ^ in the m — > limit would signal 
a spontaneous breaking of the U(4) symmetry. In the case (^a^p) oc 5 a p the symmetry 
is reduced to U(2) (8 U(2) and in the massless limit m — > there will be eight massless 
Goldstone pseudoscalar modes, parameterizing the coset U(4)/U(2) (8 U(2), with a low 
energy dynamics described by the corresponding chiral perturbation theory. The formation 
of the chiral condensate requires a sufficiently strong value of g, so there is a phase boundary 



at which the condensation "turns on." The works 
using lattice Monte Carlo methods. 



4] have located this phase boundary 
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III. DISCRETIZATION 



The fermionic part of the action ([T]) is easily discretized using the staggered fermion 
formulation [10]. The gauge field part of the action can be discretized in two ways, compact 
and noncompact, both of which will be described and used here. From this point on we 
work in units where ft, = c = 1, and use a lattice spacing at in the time direction and a s in 
the spatial directions. Thus we have lattice fields at the sites t = a t n , x-i = a s rii (i = 1, 2, 3) 
where no, • • ■ ,n% are integers. We are permitting a s ^ a t because the anisotropy parameter 
a s /a t will provide us with the handle to remove the Fermi velocity vp from the lattice 
action, so that the only parameters that will appear are the coupling ([H]) and the fermion 
mass (which must eventually be taken to zero). This mirrors the continuum redefinition 
xq = Vpt which appears in (TjJ. The lattice action takes the form 

S = \ V a t a 2 s \-{x{n)U{n)x{n + 0)-x{n)U{n-6)x{n-0)) 

noniri2 

1 2 

+v F — ^r}i{ri)(x{n)x{n + i) -j((n)x(n-i)) +mx{n)x{n) 



tE<E F^ d3) 



2 , - 

nQ...nz i=l x ' 

The notation employs four-vectors n = (n , n 1; n 2 , n 3 ) and unit vectors = (1,0,0,0) 
etc. Here x-, X are 1-component fermions and as site dependent coefficients one has the 
staggered phase factors fti(n) = (—1)™° and i] 2 (n) = (— l) n o+ n i. The reason that one- 
component fermions can be used is because staggered fermions "suffer" from doubling, 
so that in three dimensions there are eight continuum modes, which organize themselves 
into two four-component fermions under a change of basis. 1 The link fields are defined as 
U(n) = exp(iea t 6(n)), where 9(n) is the lattice version of the scalar potential A (x). Here 
we have used the noncompact form of the gauge action in the last term. The compact form 
will be discussed at a later point below. 



1 This is analogous to the four flavors that appear in the 3+1 dimensional staggered formulation of lattice 
quantum chromodynamics. 
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We next rescale to dimensionless fields, x ~~ >* x/ a s and 6 — > 9/a t e to obtain: 
S = \ E lx(n)U(n) X (n + 6)-x{n)U(n-6)x(n-6) 

n nin 2 ^ 
2 

+v F —^2r)i(n)(x(n)x(n + i) -x(n)x(n-i)) + ma^n)^) 



i=l 

E X>w-^-?)) 2 (14) 

n ...n 3 i=l 



Finally, we can absorb the Fermi speed vf into the anisotropy parameter, choosing a s /at 
vp, to obtain the lattice action in its most convenient form, 

S = \ E lx{n)U(n) X (n + b)-x(n)U(n-0)x{n-0) 

n nin 2 ^ 
2 

+ E Vi( n )(x(n)x(n + x(n)x(n - ?)) + m%(n)x(n) 



i=l 

3 



+ E f Y^(0(n) - 6(n - i)) 2 (15) 

no...n3 i=l 



where 



(j = — = — — m = ma t . (16) 
g z e z 

A slightly different choice for the anisotropy parameter a s /a t will be made below when we 
come to tadpole improvement. 

We also consider the case of a compact gauge action, where the last term in ffT5l) is 
replaced by 

3 

-0 J2 E Re u ( n ) u *( n -i) ( 17 ) 

no. ..ns t=l 

In the weak field limit (small 0(n)), which corresponds to large /3, the two formulations are 
equivalent. However, at small /3 it is expected that there will be qualitative differences. 

IV. FLAVOR SYMMETRY VIOLATION 

As stated above, a single staggered fermion automatically yields two flavors, since the 
staggered formulation does not fully solve the doubling problem. In the continuum, the 
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massless theory with two flavors has a U(4) flavor symmetry, which is reflected in a degener- 
acy of the spectrum of the Dirac operator. On the other hand, it is known that the leading 
order spectral degeneracies of the lattice Dirac operator are broken by flavor violating higher 
order terms (in the lattice spacings a t , a s ). In the massless limit but at nonzero lattice spac- 
ing only a?7(l)<8>Z7(l) flavor symmetry remains (in addition to some discrete symmetries) 



ljj . Long ago it was shown in the 3+1 dimensional case that the flavor symmetry brea king 



can be seen by going to the "flavor basis' : 



12]. For 2+1 dimensions, see for example 13] 



where staggered fermion flavor breaking terms were previously considered in the context of 
the Thirring model. Thus in the present article we are reiterating concerns that were raised 
already in 13|, though here our principal concern is the effect in the context of graphene 
effective lattice field theory. Although the flavor symmetry breaking terms are irrelevant 
operators (i.e., they are suppressed by a t ,a s ), at one loop and at finite lattice spacing they 
have important effects on the self energy of the fermions 11]. The effect of this flavor 
symmetry violation on the order parameter (ipip) that is used to distinguish the semi-metal 
versus insulator phases is not known, though in our Conclusions we will make a conjecture 
for what might occur. The flavor changing interactions are a lattice artifact that is known 
to disappear in the continuum limit. Hence, if one could send the lattice spacings a t , a s of 
the discretized effective theory (not to be confused with the lattice constant of the graphene 
system itself) to zero, one would recover the full U(4) symmetry However, the Monte 
Carlo simulations are performed at finite at,a s , and so this lattice artifact must be taken into 
account. Thus it is not quite accurate to say that one is simulating the effective theory with 
two (l+3)-dimensional Dirac fermions constrained to a plane, equivalent to four massless 
(l+2)-dimensional Dirac fermions. An extrapolation in the lattice spacing or suppression 
of the lattice artifacts is needed. One would like a systematic way to remove these lattice 
artifacts. This motivates the present study. 

We determine the size of the flavor-splitting by studying the eigenvalues of the lattice 
Dirac operator, which is the discretization of ( fl2l) corresponding to the lattice action (TI5]) . In 
Fig. [1] the "unimproved" data shows the average spectrum of the staggered Dirac operator, 
for the lowest lying modes. Here a Monte Carlo simulation was performed with /3 = 0.11, 
and eigenvalues were obtained for each configuration of the gauge field. The error bars in 
the figure indicate the standard deviation in the eigenvalues. It can be seen that there is 
a linear rise in eigenvalues, with no degeneracies whatsoever. Thus at strong coupling the 
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flavor symmetry of the continuum is badly broken. 

Next we consider the case of weak coupling, = 4.0. In Fig. [3] the unimproved data 
does show evidence of approximate degeneracies. The weaker coupling leads to smoother 
configurations of the gauge field. Rough gauge fields are farther away from the continuum 
limit, so that the 0{a t ,a s ) flavor symmetry violations is more pronounced. 

We have examined the spectrum for other values of /3. The general pattern is that for 
strong coupling the flavor symmetry is badly broken. Our next task is to attempt to restore 
it, since the j3 corresponding to graphene and the phase transition of the effective theory is 
at a strong coupling value. 



V. IMPROVEMENT 



In fact, some time ago the lattice QCD community set aside unimproved staggered 
fermions due to unwanted lattice artifacts. Modern staggered fermions are improved in 



various ways in order to suppress these effects 
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18j . So-called AsqTad staggered fermions 



were popular for several years for the study of K and B physics (e.g. 1 19]). Further improve- 
ments have been introduced to produce HISQ staggered fermions |20|. Detailed studies of 



the low lying eigenva 
been conducted in 



ue spectrum of various staggered Dirac operators have for instance 
211 ] . In each case, an important effect is to restore the flavor degeneracy 
by suppressing flavor changing interactions. The present work represents a first attempt in 
that direction; however, we will find that improvement of staggered fermions in the present 
context is more difficult. The reason is that for the study of graphene and the phase tran- 
sition of the effective theory the coupling is strong, where the flavor symmetry is badly 
broken. 

In lattice QCD it is known that flavor symmetry breaking can be ameliorated by making 
improvements to the lattice action that reduce lattice artifacts. An expansion in the lattice 
spacing a (or a t , a s in our case) and gauge coupling g allows for coefficients of various 
improvement terms to be determined in perturbation theory. However, asymptotic freedom 
should be important, since in that case it is clear how one makes these coefficients small in 
matching onto the desired continuum theory. It is then an important question whether for 
the strongly coupled theory of graphene, where there is no asymptotic freedom, the lattice 
action can be improved so as to reduce the flavor symmetry breaking effects. Certainly 
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perturbative improvement is out of the question. 



A. Tadpole improvement 

Tadpoles arise from (Aq(x)) ~ (6 2 (n)) ~ 1/af, where the estimate is made on dimensional 
grounds. As mentioned above, we study both the compact and noncompact gauge actions. 
In the noncompact case, gauge field tadpoles only enter the perturbation series through 
the gauge links U(n) = exp(ia t e9(n)) that are contained within the fermion action. In the 
compact case there are additional multiphoton vertices coming from expansion of the gauge 
action fflTj) . Consider the following example in the fermion time- like hopping terms. In this, 
we reintroduce dimensions and canonical kinetic term for 9{n) through 9{n) — > atg9(n). 
Then expanding the link U(n) = exp[ia t g9(n)] and focusing on the contribution to the 
fermion self-energy, we obtain a term a t g 2 {9 2 (n))x(n)x{n + 0) ~ {g 2 / dt)x{ n )x{ n + 0)- I- e > 
there is a large correction to the hopping term, even though the 9 2 xx vertex is irrelevant by 
power counting. There is also a large effect on the marginal 9\X vertex: 

ig9{n) (l - ^a 2 g 2 (9 2 (n)) + ■ ■ ■ J x(n)*(n + 6) 

= igOin) (1 + 0(g 2 )) x(n) X (n + 0) (18) 

Here again, the correction is 0(g 2 ) rather than 0(g 2 a 2 ), due to the tadpole (9 2 (n)) ~ 
I /a 2 . The tadpoles associated with the irrelevant fermion vertices thus give significant 
contributions to the renormalization of g, causing the matching onto continuum perturbation 
theory to be problematic. This can be circumvented through a change in renormalization 



scheme, known as tadpole improvement 15]. In fact, since for graphene the value of g 



is large, the tadpoles corrections are out of perturbative control and must be evaluated 
nonperturbatively. 

We will now show that the translation between the bare lattice = 1/g 2 (i.e., the 
parameter that appears in the action that is simulated) and its tadpole improved value /3ti 
is somewhat different depending on whether the compact or noncompact form of the gauge 
action is used, more so at stronger values of the coupling. From this, the physical coupling — 
as estimated by the tadpole improved value /3ti — is different from the bare coupling, due 
to radiative effects. In fact, we will reproduce the results of regarding the relationships 
between /3 TI and (3. 
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We begin with the expectation value (P) of the plaquette operator P = U (n)U*(n+i), i = 
1 or 2, which is related to (Aq(x)). The average link uq is defined through this quantity: 

u = (Pf 2 (19) 

Note that the square root is used here, in contrast to the fourth root that appears in 
QCD applications, since the plaquette operator is quadratic in the links that are allowed to 
fluctuate in the present, nonrelativistic formulation. 

Tadpole improvement Q can be understood as integrating out ultraviolet modes of the 
link operator U(x), to obtain an effective infrared link operator. The quantity uq represents 
the ultraviolet divergent effects of tadpoles (A^). Thus, the link is related to an infrared 
(IR) field V(n) or 6 IR {n): 

U(n) = u V(n) w u (l + ia t e6 IR (n)) (20) 
When the lattice is formulated using instead the V(n) = U(ti)/uq links, one has 
S = J2 !^(x( n )U(n)x(n + 6)-x{n)U(n-6)x(n-6)) 

n nin 2 

2 

+v F a t a s ^ Vi( n )(x(n)x(n + x(n)x(n - 2 )) + ma ta 2 s x(n)x(n) 



i=i 

+ E a t a s € jY,(6(n)-e(n-t)r (21) 

no...n3 i=l 

The redefinition of variables is now 

X=^X, 9 =—9' (22) 
a s a t e 

One finds that a s /a t = vpUo simplifies the spatial derivative term and that the result is 
equation (IT51) except that (3 and rh are replaced by 

= u o-^r = m o/3td ^ = u m T i (23) 

Note that (3^j and mxi are what would have appeared in the lattice action had we not 
included uq in the redefinition ( 1221) . Hence these are the inverse coupling and dimensionless 
mass of the tadpole improved action. By contrast, /3 and rh are the inverse coupling and 
mass that are used in the simulation after going to the redefined variables where the action 
takes its simplest form (i.e., u does not appear explicitly). Thus in the massless limit, for 
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the noncompact gauge action, the entire effect of the tadpole improvement is to rescale the 
inverse coupling according to this equation. Something similar occurs in the compact gauge 
action case. There we have in addition a factor 1/uq in front of the gauge term, 

£ 4^X>»f(» + (24) 



no— n 3 



Here then the result is 



= -P c ti (25) 

Uq 



These rescalings of /3 agree with those found recently in 



B. Naik improvement 

The Naik [lfl fermion action improvement reduces discretization errors and when the 
tadpole improvement is also performed it is given by: 



S N = a 2 s V x(n)U-[U(n)x(n + 6)-U*(n-b) X (n-6)\ 



nonin2 



+^[U(x)U(n + 6)U(n + 20) X (n + 30) 



-U*(n - 0)U*(n - 20)U*(n - 30)x(n - 30)] 
+v F a s a t ^ Vi( n )x{n)^\ci[x{n + i) ~ xO - «)] + c 2 [x{n + 32) - x(ra - 3?)] 1 

i,non\ri2 

+a 2 s a t m ^ X.(n)x(n) (26) 



n n\ri2 



Tree level improvement makes the action 0(a 2 ) accurate by setting c\ = 9/8, c 2 = —1/24. 
Next we make the redefinitions (1221) . together with setting a s /a t = vpu as before, to 
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obtain: 



S n = E x(n)Uc 1 [U(n)x'(n + 6)-U*(n-6) X '(n-6)] 
+-^|[C/(x)C/(n + 6)C/(n + 2 6)x'(ri + 3 0) 



-^-0)^- 2 0)^-3 0)^-3 0)] 

+ ^('Ox'W^ ] c i[x'( n + ? ) + c 2 [x'( n + 3? ) -x'(^-32)] > 

+m ^ x'(n)x'(n) (27) 

noniri2 

where again, m = ma t UQ. This is the "Naik-tadpole improvement;" note that Uq appears 
explicitly in this action. To obtain just the Naik improvement, one can set Uq = 1 in the 
previous expressions. 

C. Spectrum results 

We have computed the low lying eigenvalues of the spectrum of the Dirac operator on 
dynamical configurations at various values of /3, in order to see the size of the flavor symmetry 
violating effect. Fig. [1] shows the spectrum of average eigenvalues for = 0.11, rh = 0.01 
on 12 3 x 8 lattices, with compact gauge action, as well as the standard deviation (by error 
bars). Fig. [2] shows the same thing except that the noncompact gauge action was used. In 
either case, one can see that there is no hint of the four-fold degeneracy of the continuum 
theory and that the splitting is of the order 0.02. By comparison, the explicit mass in the 
simulations of [l| ranged from 0.0025 to 0.02. Thus the flavor changing interactions split the 
spectrum at the order of the mass or greater, and one is far from the desired theory. Since 



according the to Banks-Casher relation 24| the condensate on the lattice is determined by 
the density of near-zero modes, a significant systematic error will be introduced by the flavor 
splitting that we observe. We note that for the "improved" Dirac operators the splitting is 
not at all improved. This would seem to indicate that the lattice is actually quite coarse, so 
that suppressing lattice artifacts cannot be achieved by simple power-counting in the lattice 
spacings at,a s , such as is done in the Naik improvement. It is also worth mentioning that 



large scaling violations were seen in [l| for strong coupling (very small values of f3) which 
would be a further indication that lattice artifacts are playing a dominant role. However, the 
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FIG. 1: Spectrum of lowest lying modes of the three massless Dirac operators we consider, for 
compact gauge action. The configurations of gauge fields were dynamically generated at (3 = 0.11 
and m = 0.01 on a 12 3 x 8 lattice with the unimproved staggered fermion action and plaquette 
gauge term. The tadpole improvement of the Naik Dirac operators used uq = 0.256. Average 
eigenvalues are shown, and the error bars represent standard deviations. 
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» Naik 

□ Naik+tadpole 



FIG. 2: Same as Fig. [T] except that here we use the noncompact gauge action. 

n 

fact that |l[ observe scaling in a regime where we see large flavor violations is interesting, as 
it suggests that there is a universal description but that it is one with less flavor symmetry 
than the U (4) of the target graphene effective theory. 

As a further check, we have also computed the spectrum from a simulation at the weak 
coupling (3 = 4 where the flavor violation is expected to be small due to weak interactions. 
We also note that at this weak value of the coupling the compact and noncompact formu- 
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FIG. 3: Similar to Fig. [T] (compact gauge action) except that f3 = 4.0 and uq = 0.974. 

lations of the gauge action are completely equivalent. Thus the flavor symmetry breaking 
that we next describe is universal. At large ft the fluctuations in the gauge field strength 
are suppressed and a perturbative expansion of the link operators U (x) ~ 1 + iagA^ix) 
should be valid. Results for the low lying eigenvalues of the three types of Dirac operators 
are shown in Figs. [3] and HI and these certainly show a closer approximation to the four-fold 
degeneracy. However, the improved Dirac operators do not show any superiority to the 
unimproved one. This somewhat surprising result suggests that a further improvement may 
be needed, such as smearing, something which we explore in the next section. 

It is also interesting to have a statistical measure for what happens to flavor symmetry 
over an ensemble. For this purpose we have computed 



where \ is the average value of the ith. eigenvalue. This measures the relative flavor sym- 
metry breaking in the first four eigenvalues. On the /3 = 0.11 (compact) lattice we obtain 
i? FSV m 1.4(1). On the (3 = 4.0 (compact) lattice we obtain -Rfsv ~ 0.18(2). These results 
are independent of the improvement, which is curious at the larger (3. 

D. Relation between /3's 

Above, we found that the value of (3ti m the tadpole improved action can be related 
to another value /3 obtained after redefinitions, given by Eqs. ( |23|) and ( 1251) . The latter 



i(A! + A 2 + A3 + A 4 ) 



(28) 
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FIG. 4: Similar to Fig. [2] (noncompact gauge action) except that /3 = 4.0 and uq = 0.974. 



should be used in the simulation with an action that is equivalent to one without tadpole 
improvement (or only a factor 1 /uq on the temporal Naik term) . We are therefore interested 
in the effective value 0xi as a function of so that we know how to interpret simulations 
done at in terms of the underlying 0ti- For instance, Drut and Lahde find a critical value 
of the coupling for which a condensate forms, and this should be interpreted as a value of 
at which the simulation is done (i.e., in an action without Uq appearing explicitly). To see 
what this physically corresponds to, one must translate back to 0^ i n order to find the value 
of the coupling in the tadpole improved action, where ultraviolet artifacts are minimized. 

Results for the compact action are summarized in Table Hand for the noncompact action 
in Table [TTJ Thus to simulate graphene, which has 0ti ~ 0.037, we should choose the 
modified values given in the first rows of Tables [I] or El depending on the form of the 
gauge action. This gives ~ 0.12 for compact and ~ 0.004 for noncompact. The 
simulation coupling where Drut and Lahde have found a phase transition is C 0.074. 
The physical value of the inverse coupling is then approximately 0j^ 0.21, which is at 
a coupling significantly weaker than graphene, 0ti ~ 0.037. Thus the appearance of the 
condensate (ipip) occurs for a weaker value of the coupling, and will persist at the stronger 
value of graphene. One concludes that the lattice simulation is indicative of an insulator 
phase. This is in agreement with the findings of j^]. 

We also mention in passing that the value of (P) and hence uq turned out to be essentially 
independent of which fermion action (unimproved, Naik improved or Naik-tadpole improved) 
we used in the simulation. We also changed the mass to 0.02 and find the same value of uq. 
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p 


<^> 






0.037 


0.0306(40) 


0.175(11) 


0.00647(42) 


0.058 


0.034(3) 


0.183(7) 


0.01061(41) 


0.11 


0.066(4) 


0.256(8) 


0.02816(88) 


0.15 


0.0901(44) 


0.3002(73) 


0.0450(11) 


0.25 


0.1492(43) 


0.3863(56) 


0.0965(15) 


0.5 


0.504(6) 


0.710(4) 


0.355(2) 


1.0 


0.814(3) 


0.9023(14) 


0.9023(14) 


2.0 


0.9120(13) 


0.9550(7) 


1.91(14) 


4.0 


0.949(2) 


0.974(1) 


3.896(4) 



TABLE I: The average plaquette (P) and the tadpole correction factor no that is derived from it, 
as a function of f3, for the compact gauge action. This then gives a value for tadpole improved 
inverse coupling (3 = /J™. For instance, for graphene we want /3£j = 0.037 and the inverse coupling 
that should be used in the simulation is j3 0.12. 



VI. STOUT LINK SMEARING 



We have seen that at weak coupling (large /3), the spectrum degeneracies start to appear. 
This is the result of the fact that in this regime the gauge fields are smooth, whereas at 
strong coupling the gauge fields are rough. Clearly what is needed at strong coupling is a 
way to smooth out the short distance (unphysical) roughness without destroying the long 
distance (physical) fluctuations of the gauge field. The way that this can be done is to use 



smeared links in the fermion action. Here we will study stout link smearing 25| and will 
find that it successfully restores the level degeneracies for moderate to weak coupling, but 
that it fails at couplings as strong as graphene, (3ti = 0.037. 

Stout link smearing in the present context introduces the definitions 

3 

y \U(n + t) + U(n - fll. fl(n) = C(n)U*(n). O(x) = 



C(n) = pJ2[U(n + i) + U(n-i)}, fi(n) = C(n)U*(n), Q(x) = -p*(n) - fi(n)] (29) 



and U^ k '(n) at smearing step k are mapped into U^ k+1 \n) according to 



( n ) = exp [iQ( k) (n)]U {k) (n) (30) 
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<^> 


u 


^TI 


0.002 


0.0131(37) 


0.114(16) 


0.0175(25) 


0.004 


0.0108(54) 


0.104(26) 


0.0385(96) 


0.005 


0.0121(41) 


0.110(19) 


0.0455(78) 


0.01 


0.0121(27) 


0.110(12) 


0.091(10) 


0.02 


0.0118(42) 


0.108(19) 


0.184(33) 


0.037 


0.0272(46) 


0.165(14) 


0.224(19) 


0.058 


0.0757(40) 


0.2751(73) 


0.2108(56) 


0.11 


0.2392(45) 


0.4891(46) 


0.2249(21) 


0.25 


0.5228(40) 


0.7230(27) 


0.3458(13) 


0.5 


0.7192(29) 


0.8481(17) 


0.5896(12) 


1.0 


0.8466(20) 


0.9201(11) 


1.0868(13) 


2.0 


0.9195(11) 


0.9589(6) 


2.0857(13) 


4.0 


0.94(2) 


0.97(1) 


4.124(43) 



TABLE II: The average plaquette (P) and the tadpole correction factor uq that is derived from it, 
as a function of f3, for the noncompact gauge action. This then gives a value for tadpole improved 
inverse coupling f3 = /J™. For instance, for graphene /S™ = 0.037 and the inverse coupling that 
should be used in the simulation is /3 ~ 0.004. 

It can be seen in Fig. [5] that smearing works very well at weak coupling. The smeared 
eigenvalue data has 10 smearing iterations with smearing parameter p = 1/6, where the 
latter was found to be optimal based on trial and error. Less smearing iterations obviously 
results in less degeneracy. Unfortunately, as the coupling is made stronger, the smearing 
becomes progressively less effective, as can be seen in Fig. [TJ 

By contrast for the noncompact gauge action, even at the relatively small value of 
(3 = 0.11, one finds a significant improvement from smearing; see Fig. |6j Since the phase 
transition occurs at j5 « 0.07 we expect smearing to be quite useful for reducing flavor sym- 
metry breaking in the vicinity of this point. On the other hand from Table HT1 we found that 
graphene with /3jPj = 0.037 corresponds to /3 « 0.004 which is far too strong for smearing to 
help. Indeed we have found that there is no restoration of degeneracy in this case. 
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FIG. 5: P = 4.0, compact gauge action, 10 stout link smearings with p = 1/6. 
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FIG. 6: P = 0.11, noncompact gauge action, 10 stout link smearings with p = 1/6. 



VII. CONCLUSIONS 



We have found that at P ^ 1 both the unimproved action, and the tadpole-improved 
action with a Naik term show significant flavor symmetry breaking. We have also measured 
the average plaquette term used for tadpole improvement and have described how it calls for 
a reinterpretation of previous lattice simulation results. Importantly, it indicates that the 
insulator /semi- metal phase transition observed on the lattice occurs at a physical coupling 
that is significantly weaker to the one that appears in suspended graphene. It follows that 
the lattice simulations predict that the chiral symmetry is spontaneously broken and that 
suspended graphene would be in the insulating phase. 

On the other hand, conjectured phase diagrams in the g versus Nf plane would indicate 
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that the critical g decreases as Nj is decreased. So, if the staggered formulation really 
simulates effectively less that Nf = 2 due to the flavor symmetry breaking, the lattice 
simulations would predict a critical g that is weaker than that of graphene. Restoration of 
the U(2Nf) flavor symmetry would tend to increase the value of the critical g. Thus it is still 
possible, though unlikely, that lattice simulations would predict that suspended graphene 
is in the semi-metal phase, provided the full flavor symmetry is intact. We think that it is 
unlikely since the critical ft c would have to shift all the way from ft = 0.07 to ft — 0.004 as 
a result of restoring the flavor symmetries. Still, a study with overlap fermions is of interest 
to settle the question. 

We have conducted studies with both the compact and noncompact formulation in their 
gauge action. In 1+3 dimensional quantum electrodynamics, the compact formulation has 
difficulties with a bulk phase transition in the strong coupling regime, separating it from 
the continuum theory (see for example 22] and recent work in 23]). On the other hand, 
with the nonrelativistic constraint Ui(x) = 1, i = 1,2,3 that we impose, the phase structure 
of the compact theory will be quite different since, for instance, magnetic monopoles will 
not exist. However, the presence of vortices requires further investigations of the compact 
theory, which we will leave to future work. At present what is known from [3!] is that the 
compact theory has a first order phase transition in contrast to the second order transition 
of the noncompact case. This very different phase structure indicates that nonperturbative 
features, such as vortices, are having a significant effect in distinguishing the two theories at 
strong coupling. In the present article we show results for both compact and noncompact 
gauge action. We find that the qualitative features do not change: the large flavor violations 
are present in either formulation at strong coupling. 
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Appendix A: Simulation details 

All of our results were obtained using hybrid Monte Carlo simulations with dynamical 
staggered fermions. This simulation method has been reviewed in the present context in 
Q]. The mass in our simulations was ma = 0.01, where a is the lattice spacing. We have 
simulated on various sizes of lattices (6 3 x 8, 8 3 x 8, 12 3 x 8, 16 3 x 8 and 24 3 x 8). We checked 
that the configurations were fully thermalized by comparing ordered and disordered starts. 
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